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Scaling Relations from Sunyaev-Zel’dovich Effect and Chandra 
v X-ray measurements of high- redshift galaxy clusters 

Massimiliano Bonamente 1 ’ 2 , Marshall Joy 2 , Samuel J. LaRoque 3 , John E. Carlstrom 3,4 , 

Daisuke Nagai 5 and Dan Marrone 6,3 

ABSTRACT 

We present Sunyaev- Zcl' dovich Effect (SZE) scaling relations for 38 massive 
galaxy clusters at redshifts 0.14 < z < 0.89, observed with both the Chandra 
X-ray Observatory and the centimeter- wave SZE imaging system at the BIMA 
and OVRO interferometric arrays. An isothermal 3-model with central 100 kpc 
excluded from the X-ray data is used to model the intracluster medium and 
to measure global cluster properties. For each cluster, we measure the X-ray 
spectroscopic temperature, SZE gas mass, total mass and integrated Compton-y 
parameters within r 25 oo- Our measurements are in agreement with the expec- 
tations based on a simple self-similar model of cluster formation and evolution. 

We compare the cluster properties derived from our SZE observations with and 
without Chandra spatial and spectral information and find them to be in good 
agreement. We compare our results with cosmological numerical simulations, and 
find that simulations that include radiative cooling, star formation and feedback 
match well both the slope and normalization of our SZE scaling relations. 

1. Introduction 

The Sunyaev-Zel‘ dovich Effect (SZE) is a unique and powerful observational tool for 
cosmology (for review see Carlstrom et al. 2002). It is a small distortion in the cosmic 
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microwave background (CMB) spectrum caused by scattering of CMB photons off a distri- 
bution of high energy electrons in dense structures such as clusters of galaxies (Sunyaev k 
Zel’dovich 1970, 1972). This effect has a unique property that the signal is independent of 
redshift, making it particularly well suited for deep cluster surveys (e.g., Holder et al. 2000; 
Weller et al. 2002). Several Sunyaev-Zel'dovich Effect survey experiments are currently in 
progress (Ruhl et al. 2004; Fowler 2004; Kaneko 2006), and are expected to generate a large 
sample of SZ-selected clusters with masses greater than ~ 2 x 10 w M s . The resulting large 
samples of galaxy clusters will enable direct measurements of the evolution of the number 
density of galaxy clusters as a function of redshift and in principle can provide a powerful 
constraint on the nature of dark energy (Wang k Steinhardt 1998; Viana k Liddle 1999; 
Mohr et al. 2000; Haiman et al. 2001). 

To utilize the upcoming SZE cluster surveys for cosmological studies, it is important 
to understand the relation between the SZE observables and the mass of the cluster. If 
the evolution of clusters is dominated by gravitational processes, a simple model of cluster 
formation and evolution based on the virial theorem (Kaiser 1986) predicts simple power- 
law relations between cluster masses and certain integrated cluster properties, including the 
integrated SZE flux (which is proportional to Y, the integral of the Compton- y parameter 
over the solid angle of the cluster). Numerical simulations further suggest that Y should be 
an excellent proxy of cluster mass when measured on sufficiently large scales (e.g., da Silva 
et al. 2004; Motl et al. 2005; Nagai 2006). These simulations also predict that the slope 
and redshift evolution of the SZE scaling relations are relatively insensitive to the details of 
cluster physics, although numerical simulations show that the input cluster physics affects 
the normalization of the SZE scaling relations (Nagai 2006). It is therefore important to 
investigate the properties of the SZE scaling relations observationally. 

Previous studies have addressed the correlation between the SZE signal and X-ray prop- 
erties. For instance, Cooray (1999) find a positive correlation between the central SZE 
decrement and the X-ray luminosity in a sample of 14 clusters. Similarly, McCarthy et al. 
(2003) detected correlations between the central SZE decrement and X-ray determined mass, 
temperature, and luminosity for a 22 cluster sample. These studies use data from multiple 
SZE and X-ray experiments, making systematics more difficult to control, and focus on the 
relationship between the central values of the SZE signal with the X-ray properties. More 
recently, Benson et al. (2004) showed that the integrated SZE flux is a more robust observ- 
able than the central values of the SZE signal, and found a strong correlation with X-ray 
temperatures using a sample of 15 clusters obtained by the Sunyaev-Zeldovich Imaging Ex- 
periment (SuZIE, Holzapfel et al. 1997; Benson et al. 2003) and X-ray temperatures from 
the ASCA experiment. 
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In this paper, the third in a series of papers combining Sunyaev-Zel’dovich Effect and 
Chandra X-ray measurements of galaxy clusters to study cosmological properties (Bonamente 
et al. 2006; LaRoque et al. 2006), we present observational studies of Sunyaev-Zel’dovich 
Effect scaling relations for clusters of galaxies. This paper advances the results of previous 
cluster scaling relations works in several ways. First, we use the largest observational sample 
yet constructed (38 clusters at redshift z=0. 14-0.89). Second, our analysis is based on SZE 
and X-ray observations obtained using the same instruments: all SZE fluxes are determined 
using centimeter- wave interferometric data from the BIMA/OVRO SZE imaging experiments 
(e.g. LaRoque et al. 2003), and all cluster X-ray properties are derived using data from the 
Chandra X-ray Observatory. Finally, our Chandra observations have an order of magnitude 
better spatial resolution than the X-ray data used in previous studies. These observations 
therefore result in more accurate measurements of the global cluster properties. 

Throughout the paper, we assume a ACDM cosmology with Qm= 0.3, Oa= 0.7 and 
h = 0.7, where h is defined such that H 0 .= 100 h km s" 1 Mpc -1 . All uncertainties are at 
the 68.3% confidence level. 

2. Theory of cluster scaling relations 
2.1. The virial radius and r 2 soo 

In order to establish relationships between mass, SZE flux and other cluster properties, 
one needs to define an outer radius - out to which all quantities will be calculated - that is 
physically motivated, reachable with the current X-ray and SZE observations, and equivalent 
for clusters of different redshift. One candidate is the virial radius. In a Friedman-Robertson- 
Walker universe, an unperturbed spherical region expands indefinitely, while a perturbed 
overdense region (the seed of a future cluster) eventually recollapses. When the overdense 
region collapses under the effect of its own gravity, it is assumed to reach virial equilibrium 
when the radius is half of that at maximum expansion (Peebles 1980; Lacey & Cole 1993). 
The ratio of the mean cluster density to the background density at the time of virialization 
is A v = 187r 2 for a universe with critical matter density (0l m = 1). For a different cosmology 
with fife = 0, Bryan & Norman (1998) showed that A r ~ 187T 2 + %2x — 39a; 2 , where x - 
fl m o(l + z) 3 /E 2 (z) and E 2 {z) = fl m o(l + z ) 3 + fh\ + f2*o(l + z) 2 , as found from a fit to 
numerical simulations (Lacey & Cole 1993). 

With this characterization of the mean cluster density at time of virialization, the virial 
radius can be determined as the radius within which the average density of the cluster is A„. 
times the critical density, via 
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-■Kp c {z)A v (z)rl ir = M tot (r^ {z) ) 


( 1 ) 


in which both p c (z) and A v (z) are cosmology dependent, and the critical density p c {z) is 
defined as: 

( 2 ) 


Unfortunately, the virial radius is usually unreachable with current X-ray and SZE 
measurements, and one is forced to perform measurements out to a smaller radius. Such 
radius (r a) is characterized by the density contrast parameter A in place of A^z) in Equation 
2, and corresponds to a higher average density, 4/3 • 7rp c (z)A • r\ = M tot {r/\). We choose a 
contrast parameter A = 2500, corresponding to an average density of 2500 times the critical 
density at the Cluster’s redshift. This choice is motivated by the fact that this is the radius 
that is typically reachable with our Sunyaev-ZePdovich Effect and X-ray data without any 
extrapolation of the models (Bonamente et al. 2006; LaRoque et al. 2006). 7 


2.2. Scaling relations 

The hierarchical structure formation theory developed by Kaiser (1986) predicts simple 
relationships between physical parameters of collapsed structures, known as scaling relations. 
With the assumptions of hydrostatic equilibrium and of an isothermal distribution for both 
the dark matter and the cluster gas (e.g., Bryan & Norman 1998), it can be shown that there 
is a simple relationship between a cluster’s total mass and its gas temperature T e : 

T e oc M% t 3 E(z) 2 ' 3 (3) 

where the mass was calculated out to a radius of mean overdensity A, M to t = M tot (r a). 
For f gas = M gas /M tot , the expected relationship between the gas mass within ?’a ( M gas = 
Mgas(r a)) and the gas temperature is 

T e ff a l a M% 3 s E(z) 2 / 3 . (4) 


7 The use of a constant overdensity factor A was shown by Maughan et al. (2006) to give results similar 
to the case of a variable overdensity factor A(z) — A(0)[A v (^)/A v (0)], in which the variable overdensity 
scales with redshift in order to keep the ratio of two comoving densities constant. 
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The Compton-y parameter is a measure of the pressure integrated along the line of sight: 

y = / &T n e o dl ( 5 ) 

J 0 "TeC 2 

One can further integrate the y parameter over the solid angle Q subtended by the cluster, 
to obtain the integrated Compton-y parameter: 

<6) 

where A is the area of the cluster in the plane of the sky. In the context of an isothermal 
model, Y is proportional to the integral of the electron density n e over a cylindrical volume, 
thus 

YD 2 a (x T e J n e dV = M gas T e = f gas M tot T e . (7) 

In section 5 we consider the effect of integrating gas mass within a spherical volume while 
determining Y in a cylinder. Using equation 3 we can rewrite equation 7 in terms of either 
M tot or T e , or substitute M gas /f gas for M toU to obtain: 

f 9 a S T? l2 E(z) ' 

( 8 ) 

JgaT M glsE(z)V 3 

Equations 8 are the scaling relations that we investigate observationally in this paper. 

3. Sunyaev-ZePdovich Effect and Chandra X-ray observations of galaxy clusters 

3.1. Data 

We analyze the Sunyaev-ZePdovich Effect and X-ray data observations of 38 clusters in 
the redshift range z=0. 14-0.89, observed with the BIMA and OVRO interferometric arrays, 
and with the Chandra X-ray imaging spectrometers. The data are presented in Bonamente 
et al. (2006) and LaRoque et al. (2006), the previous two papers in this series, in which the 
observations were used respectively to measure cluster distances and the Hubble constant, 
and the X-ray and SZE gas fraction. We refer to LaRoque et al. (2006) for details on 
the observations and data modeling. In the following, we review those aspects of the data 
modeling and data analysis that are relevant to the investigation of the scaling relations. • 
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3.2. Data modeling 


The gas density model is based on the spherical /3-model (Cavaliere & Fusco-Femiano 
1976, 1978), which has the form 

/ r 2\ -3/J/2 

n e(r) = n e o ^1 T 


where n e o is the central electron number density, r is the radius from the center of the cluster, 
r c is a core radius, and fi is a power-law index. When integrated along the line of sight to 
determine the projected SZE decrement distribution (oc n e ), and X-ray surface brightness 
(oc rig) this model has the simple analytic forms 

I AT 
S* 

s, 

where A T 0 is. the central thermodynamic SZE temperature decrement /increment and 9 C 
is the angular core radius of the cluster (e.g., Birkinshaw et al. 1991; Reese et al. 2002). 
This model typically provides a good description of the X-ray surface brightness and SZE 
decrement profiles out to ~r 250 o (e.g., Jones & Forman 1984; Elbaz et al. 1995; Grego et al. 
2001; Reese et al. 2002; Ettori et al. 2004). This simple model, however, does not provide 
a good description of the peaked X-ray surface brightness observed in the center of some 
clusters. To minimize the systematic bi;us associated with modeling of cluster cores, we 
therefore exclude the central 100 kpc from both the spatial and spectral X-ray data, as was 
done in the previous two papers in this series. The X-ray spectroscopic temperature is also 
determined by a single-temperature (it to the X-ray spectrum of photons extracted from an 
annulus between 100 kpc and r 2 5 oo (LaRoque et al. 2006). This 100 kpc-cut model was shown 
to recover the gas masses of simulated clusters with a range of dynamical states to better 
than 5% accuracy at r 2500 (LaRoque et al. 2006) and, when applied to the determination 
of the Hubble constant, yielded the same results as a more complex non-isothermal model 
(Bonamente et al. 2006). 


AT 0 \1 + 

S x 0 ( 1 + XX ) 


02 V (l-3/3)/2 

W 

0 2 \ (l-6/3)/2 

WJ ’ 


(9) 


3.3. Analysis Methods 

Best-fit model parameters and confidence intervals for all model parameters are obtained 
using a Markov chain Monte Carlo (MCMC) method described in detail by Bonamente et al. 
(2004, 2006). LaRoque et al. (2006) explains the implementation of the likelihood calculation 
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for the 100 kpc-cut model. For each cluster, the Markov chain constrains the parameters 
5*0, 0 ci A T 0 , T e , and A. We use the cosmological parameters h = 0.7, Qm— 0.3 and Q,\=0.7 
to calculate each cluster’s angular diameter distance Da (e g., Carroll et al. 1992). 

From these model parameters we calculate r 25 oo and A7/„<(/250o) through the hydrostatic 
equilibrium equation (e.g., Grego et al. 2001), 


M tot (r ) 


3 pkTx r 3 
Gpm p r% + r 2 ’ 


( 10 ) 


in which fi is the mean molecular weight calculated using the X-ray metallicities and r c = 
Q c Da ■ One obtains r 2 5 oo from the solution of the following equation: 


wm rl 
G/mip r% +r\ 


/l q . . . 

Y r A A Pc(z), 


in which p c (z ) is given by Equation 2, A = 2500, and the right hand side is just Mt ot (r 25 00 ). 
We then compute the global cluster quantities needed for the analysis of SZE and X-ray 
scaling relations. The gas mass is computed by integrating the gas density model, 


M gas (r 25 oo) = AKp e n e0 m p D A 


r a ( 


rp 
1 + K 


-3(9/2 


e 2 de, 


(ii) 


where p e is the mean molecular weight of the electrons, and n e0 is the central electron density, 
obtained from the parameters of the (3 model (LaRoque et al. 2006, Equation 12). 

The integrated y parameter (Y, equation 6) is calculated using the measured SZE decre- 
ment AT, which is directly proportional to the Compton-y parameter 


AT = Tcmb ■ f {x) ■ y. 

The factor f{x) is the frequency dependence of the SZE: 

f{x) = (x ■ coth (!) - 4^ • (1 +< S W), 

in which a; = hv/kBTcMB, and 5 re i is a small relativistic correction factor. At our observing 
frequencies, f(x) ~ - 2. Thus, 


Y = 


L 


AT 


TcmbJ (x) 


dQ = 


at, 


rr 2 50o/D A / q2 \ 

Jo ( l + ?) 


l CMB 


/(*) 


( 1 — 3 / 3)/2 


Odd. 


( 12 ) 
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4. Observational constraints on Sunyaev-Zel’dovich Effect scaling relations 

4.1. Regression method 

Our measurements of masses and integrated y parameters, using the method described 
in Section 3, are shown in Table 1. The errors in Table 1 represent the photon-counting 
statistical uncertainties of the X-ray data and the statistical uncertainties of the SZE obser- 
vations. Additional sources of uncertainty in the measurements of masses and temperatures 
include cluster asphericity, small-scale clumping of the gas, presence of point sources in the 
field, CMB anisotropy etc., as discussed in Bonamente et al. (2006) and La Roque et al. 
(2006). In fitting Y versus Tx, M gas and M tot (Equations 8), we include an additional source 
of error in the amount of ±20% for the masses and Y, and of ±10% for the temperature, as 
found in the uncertainty analysis described in Bonamente et al. (2006) and LaRoque et al. 
(2006). 

We perform a linear least-squares regression in log space, log(Y) — A ± Blog(X), 
following the method of Press et al. (1992) and Benson et al. (2004). This method accounts 
for errors in both measured parameters for each scaling relation, and it minimizes the y 2 
statistic defined as 

2 . (M^) - A ~ BlogjJQf 

a L/(Vi) I {Baiog(Xi )) 2 

in which cr loy{Yi) = a Yi /Yilog(e), (Ji og { Xi ) — oxJX, fog (e), and the linear errors a Yi and a Xi 
are obtained from the upper and lower uncertainties around the best-fit values as a = 
(<r + ± a~)/2. 

4.2. The Y -M gas , Y-M to t and Y-kT scaling relations 

The above derivation of the self-similar scaling relations does not include any variation of 
the gas fraction with cluster mass. However, there may be some evidence for such variation 
in both X-ray observations (e.g., Vikhlinin et al. 2006) and simulations (Kravtsov et al. 
2005). We examine this in the present data by performing a logarithmic fit to the f ga s~M gas 
data using a linear relationship (log{Y) = A ± Blog(X)). We find no significant evidence 
for a variation of f gas with mass (B = 0.14 ± 0.08). In the following we therefore assume 
that f gas is a constant. We then perform similar logarithmic fits to the Y-M gas , Y-Mt ot and 
Y-kT data. The results are shown in Figure d and in Table 2. Under the assumption of a 
constant f gas , all scaling relations are consistent (within 2cr statistical uncertainty) with the 
simple self-similar model of cluster evolution. 
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4.3. Redshift evolution of the Y-M gas , Y-M tot and Y-kT scaling relations 

The large number of clusters (38) and redshift coverage of our Chandra and 0VR0- 
BIMA data (0.14 < z < 0.89) enables the investigation of a possible redshift evolution of 
the SZE scaling relations. For this purpose, we divide our sample evenly into low-redshift 
clusters (z < 0.30. 19 clusters) and high-redshift clusters (0.30 < 2 < 0.89. 19 clusters), and 
repeat the logarithmic fits of Sections 4.2. 

The results of Table 2 indicate that there is no evidence for redshift evolution of the 
SZE scaling relations. First, f gas is consistent with a constant for both low and high-redshift 
clusters. Moreover, the SZE scaling relations are consistent with the self-similar slopes at 
(or below) the 2.5cr level. 

5. Comparison of Sunyaev-Zel’dovich Effect and X-ray measurements 

In the previous section we have examined the Sunyaev-Zel’dovich Effect scaling relations 
based on quantities derived jointly from SZE and X-ray observations. Here we compare 
the cluster properties derived from our SZE observations without considering the X-ray 
data, in order to determine whether the relations we observe depend strongly on the X-ray 
information. 

We analyze the SZE data using the model of Section 3.2. using additional assumptions 
to provide the constraints that would otherwise be provided by the X-ray data. As a first 
assumption we fix f3 = 0.7 and fit for 9 C and ATo in Equation 9, following LaRoque et al. 
(2006). The data quality allows us to perform this SZ-only analysis for 25 of the clusters in 
the full sample, as shown in Table 3. Knowledge of the gas temperature is required in order 
to determine r 25 oo> as can be seen from the combination of Equations 1 and 10: 



In the absence of complementary X-ray spectroscopic data, we estimate the gas temperature 
directly from the SZE data following the iterative method described by Joy et al. (2001). 
We first choose an initial estimate of the gas temperature, from which we obtain r 25 oo and 
A/t 0 t (r 2 5 oo)- We derive M gas (r 2500 ) using the equations described in Section 3.3, with the 
central gas density n e o calculated from the parameters of the SZE decrement model (Equation 
13 of LaRoque et al. 2006). We provide a final constraint by assuming that the gas mass 
fraction of each cluster is equal to the average value for this sample, f ga s = 0.116 (LaRoque 
et al. 2006), and iteratively solve the equation M ffas (r 25 oo) = f gas • M tot (r 250 0 ) in order to find 
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self-consistent estimates of T and r 2 5 oo- 

Results of the SZE analysis are shown in Table 3. In Figure 3 we compare the SZE 
measurements of gas temperature, 725005 Y, and gas mass with the values from the joint 
analysis of Section 3. The quantities inferred from the SZE data are in good agreement 
with those from the joint analysis, indicating that we have not altered the scaling relations 
of Section 4.2 by incorporating the X-ray data. Although these results show that it is 
possible to estimate Y arid the cluster gas mass from the SZE data alone, the joint analysis 
is preferred because it provides the strongest constraints on cluster properties, requires fewer 
assumptions about the cluster structure and composition, and can be applied to the full 38- 
cluster sample. Nevertheless, the ability to derive cluster properties from SZE observations 
alone will be important for the many Sunyaev-ZePdovich Effect cluster surveys that are 
currently underway (Ruhl et al. 2004; Fowler 2004; Karieko 2006), which are expected to 
generate large samples of SZ-seleeted clusters but will not generally have access to deep 
X-ray observations for cluster characterization. 


Recent work by Kravtsov et al. (2006) indicate that the quantity Y x ; = M gas * fcjgT e , the 
X-ray analogue of Y, is a low-scatter proxy for the cluster total mass. For the isothermal 0- 
model used in this paper, the gas mass can be estimated from X-ray data by using Equation 
11, and thus our data provide also a measurement of Yx. We can therefore compare Y with 
Yx, in order to establish observationally that the two quantities are indeed equivalent. In 
the case of the isothermal /Amodel, the integrated Compton-?/ parameter is an integral of 
the electron density over a cylinder G of infinite length along the line of sight, and of area 
A = rrrf 500 : 


YD\= f n e (r)dV = — f n e m pt x e dV. 

\ rn e cr JJc \ ™, e c 2 - ) m p (j, e J c 

Since the gas mass is given by an integral over a sphere S of radius r 25 oo, 



n e m v \x e diV 


the relationship between Y and Yx is 


YD\ 



(13) 


where the constant C = f c n e dV j f s n e dV accounts for the different domain of integration 
of Y and Y x , and depends on the parameters of the (3 model. We calculate C separately for 
each cluster, and typically find C ~2. 


In Figure 4 we plot Y, as derived from the joint SZE/X-ray analysis, against Y X - The 
agreement with the relationship in Equation 13 shows that Y X is indeed an unbiased esti- 
mator of the integrated Compton-y parameter. 
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6. Comparison with theoretical simulations 

We compare our results with those of recent cosmological cluster simulations (Nagai 
2006; Nagai et al. 2007), which include radiative cooling and UV heating, star formation and 
stellar feedback processes in addition to the standard gas dynamics. In Figure 5 we compare 
16 clusters simulated at z=0 and 0.6 with cooling and star formation feedback processes 
(in red), the same sets of clusters performed using non-radiative gas dynamics (in green), 
and our 38 clusters observed with Chandra and OVRO-BIMA (in black). The comparison 
indicates that both simulation models show the same slope as the observed clusters, while 
the cooling and star formation feedback model is a closer match to the observed clusters for 
the normalization. 


7. Discussion and conclusions 

Our analysis of 38 clusters observed with Chandra and OVRO-BIMA indicate that the 
SZE observable Y (the integrated Compton-y parameter) of galaxy clusters follow simple 
power-law scaling relations that are consistent with a self-similar model of gravitational 
cluster evolution. We have investigated three scaling relations between Y and total mass, 
gas mass and gas temperature. Fits of the Y-M gas , Y-M tot and Y-kT data to a power-law 
model agree with the slope predicted by a self-similar model in which the evolution of clusters 
is dominated by gravitational processes. The normalization of the Y-M gas scaling relation 
agrees well with the numerical simulations of Nagai et al. (2007), in which collisionless 
dyanmics of dark matter and gas dynamics are complemented by cooling, star formation and 
feedback phenomena. The agreement provides observational evidence that non-gravitational 
phenomena may also be an important factor in the physics of clusters. 

The redshift coverage of our sample enabled an analysis of the scaling relations as 
function of redshift, by defining allow redshift sample (0.14 < z < 0.30, 19 clusters) and a 
high-redshift sample (0.30 < z < 0.89, 19 clusters). Both samples follow the prediction based 
on the self-similar model. Our data indicate no significant evolution in the SZE properties 
of clusters at redshift z < 1 . 

We also measure the cluster mass and integrated Y parameter using the SZE data alone, 
without making use of the Chandra spectral and spatial information. These measurements 
are in good agreement with those based on the joint X-ray/SZE analysis, providing quan- 
titative evidence that SZE surveys can be used to determine the number density of galaxy 
clusters as function of mass and redshift. 
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Table 1. Cluster Parameters from Joint Analysis of X-ray and SZE Data 


Cluster z Da E(z) r 2 5oo kT M gas M tot Y f ga $ 

(Gpc) (”) (keV) (10 13 M q ) (10 14 Mq) (10- 10 ) 

Abell 1413 0.14 0.52 1.07 206±f 7.5±g-| 2.6±gg 2.2±g? 2.99±g|| 0.12±ggg 

Abell 1689 0.18 0.63 1.09 219±| 10.5±g-| 5.1±g| 5.0±g| 3.79±gff 0.10±ggg 

Abell 1835. 0.25 0.81 1.13 172±t 11.4 ±g£ 5.8±g| 5.6±g| 2.09 ±g}| 0.10±g 08 

Abell 1914 0.17 0.60 1.09 228±| 9.5±g | 4.8+gj 4.8±g| 3.01±gfg 0.10±g gg 

Abell 1995 0.32 0.96 1.18 133±| 8.2±g| 3.5±gJ 4.7±g-| 0.75±gg| 0.07±g;gg 

Abell 2111 0.23 0.76 1.12 141±£° 8.2±Jg 2.2±gl 2.5±gg 0.95±g|J 0.09±gg} 

Abell 2163 0.20 0.68 1.10 206±| 14.8±g| 8.1±g§ 5.5±g| 6.89±g gg 0.15±ggg 

Abell 2204 0.15 0.54 1.08 256+J? 11.2d=g f 4.7±g?> 5.0±gg 4.43±g|g 0.10+gg} 

Abell 2218 0.18 0.63 1.09 191±g 7.8±gt 3.0±gg 3.3±gf 1.94+gJg 0.09±g gg 

Abell 2259 0.16 0.57 1.08 172±| 5.8±g | 1.8±g} l,8±g | 0.82±g|g 0.10±g gJ 

Abell 2261 0.22 0.73 1.12 148±| 7.4±gg 3.0±g| 2.6+gf 1.34±g 0.12d=g gi 

Abell 267 0.23 0.76 1.12 132±? 5.9±gi 2.2±gl 2.0±g| 0.72±gjg O.llig gJ 

Abell 370 0.38 1.07 1.22 97±| 8.7±gg 2.8±g| 2.8±gf 0.71±g g| 0.10±ggg 

Abell 586 0.17 0.60 1.09 182±f 6.4±g | 2.3±g} 2.5±g| 1.03±g-^ 0.09±gg} 

Abell 611.... 0.29 0.90 1.16 111+f 6.8±g| 2.4±g} 2.1±gl 0.54±ggg 0.11+ggj; 

Abell 665 0.18 0.63 1.09 162±f 8.4±g-| 2.6±gj 2.0±g-} 2.68±g^ 0.13±ggg 

Abell 68.... 0.26 0.83 1.14 153±J° 9.6±{g 3.6±g| 4.3±g? 1.01±g}| 0.08+g-gJ 

Abell 697 0.28 0.88 1.15 134+f 10.2±g| 4.4±g g 3.5±g| 1.67+ggg 0.13+ggj 1 

Abell 773... 0.22 0.73 1.12 148±g 8.2 ! gg 2.7±g'l 2.6±gg 1.68±g}| 0.11±ggJ 

CL J0016+1609. .... 0.54 1.31 1.34 80±| 10.5±gf 4.4+g-f 3.3±g | 0.73±g gg 0.13+ggJ 

CL J1226+3332 0.89 1.60 1.65 66±| 13.5±|| 3.9±g| 5.2±f g 0.35±ggf 0.08±g$ 

MACS J0647.7+7015 0.58 1.36 1.37 92±g 14.1±fg 4.9±g| 6.0+}? 0.62=fcg gf 0.08+ggj; 

MACS J0744.8+3927 0.69 1.47 1.47 59±| 8.1±g f 3.1d=g | 2.3±g : | 0.34±g : g| 0.14+ggJ 

MACS J 1149.5+2223 0.54 1.31 1.34 71±| 9.9±g f 3.1±g;| 2.3±gf 0.58±gg? 0.13±gg} 

MACS J1311.0-0310. 0.49 1.25 1.30 74±f 7.2±} : ? 2.1±g : | 2.2±g : £ 0.28±g : gg O.lOzbg g| 

MACS J1423.8+2404 0.55 1.32 1.35 66±| 7.0±g| 2.3±g j; 1.9±g : | 0.28±g : g| 0.12+ggJ 

MACS J2129.4-0741 . 0.57 1.35 1.36 73±| 8.6±g g 3.3+gf 2.8±g ; g 0.39±g : g| 0.12±g gJ 

MACS J2214.9-1359. 0.48 1.23 1.29 91±| 10.2±jg 3.9+gg 3i9±g : g 0.77±g : gg 0.10+ggJ 

MACS +2228:5+2036 0.41 1.12 1.24 81±| 8.4±g | 2.8d=g‘| 2.0 hg g 0.94±gg 0.14+ggi 

MS 0451.6-0305 0.55 1.32 1.35 82±| 9.9±g f 4.8±g j 3.8±g g 0.66±g : gf 0.13±gg} 

MS 1054.5-0321 .... . 0.83 1.57 1.59 89±f 9.8+JJ 7.4±}g 4.5±J : g 0.77+ggJ 0.16±g : gi 

MS 1137.5+6625.... 0.78 1.54 1.55 42±f 4.5±gf 1.2±gg 1.0+g : | 0.09±g : gf 0.12±g : g} 

MS 1358.4+6245 ... . 0.33 0.98 1.19 113±g 8.9±g f 2.5+g-f 3.1±g ; g 0.56±g ; gf 0.08+g : g| 

MS 2053.7-0449 0.58 1.36 1.37 54 fcg 4.8±g| 0.9±g : { 1.2±g : | O.OO.t.g gg 0.08 l-g : g| 

KX .11317.5-1 145.... 0.45 1.19 1.27 122±| 16.5+Jg 8.8±g : | 8.1±g ; | 1.62±ggf 0.11±g : gg 

HX .11716.44 6708... 0.81 1.56 1.57 45±t 6.6±J j 1.2±g| 1.4+ gg 0.10±gg| 0.09 : g : g! 

HX .12129.71 0005... 0.24 0.78 1.13 128.1 g 6.7lgg 2.6.1 g ; 'f 2.1±g : | 0.66±g : fg 0. 12J.gg| 

ZW 3146 0.29 0.90 1.16 132 1 1 8.3 i g i* 1.4 Lg { 3.6±g ; | 0.88 1 g | 0.12±g gg 
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Table 2. Scaling Relations from*. Joint Analysis of X-ray and SZE Data 


Scaling 

relation 

All clusters 

0.14<z<0.30 

0.31<z<0.89 


A 

B 

A 

B 

A 

B 

Jgas i Mg (l s 

-2.86+1.09 

0.14+0.08 

-2.60+1.79 

0.12+0.13 

-3.00+1.37 

0.15+0.10 

Y,kT 

-6.17+0.23 

2.30+0.25 

-6.29+0.36 

2.42+0.37 

-6.04+0.31 

2.18+0.32 

Y,Mgas 

-22.64+1.90 

1.37+0.14 

-24.89+3.61 

1.53+0.26 

-21.10+2.21 

1.25+0.16 

Y ,M tot 

-26.39+2.90 

1.53+0.20 

-28.95+5.14 

1 .71+0.34 

-24.10+3.38 

1.38+0.23 


Table 3. Cluster Parameters from Analysis of SZE Data 


Cluster z r 2 5 oo kT M gas M tot Y 

(”) (keV) (10 13 Mq) (10 14 Mq) (10- 10 ) 

Abell 1689 0.18 196±f 8.0±gf 4-2+gf 3.6+g\f 1.88±gfg 

Abell 1835.......... 0.25 169±| 11.4+jg 6.1±g : f 5.3±g ; | 2.6(H g’fg 

Abell 1914.......... 0.17 212±g° 8.5+i| 4.5±g : | 3.9- Iff; 2.39±g : n 

Abell 1995...... .... 0.32 117±f 8.4±g : g 3.7±g : f 3.2±g| 0.84±g ; Jf 

Abell 2111. . 0.23 124±g 5.4±j : g 2.0±g : f 1.7±g : | 0.45±g ; ?l 

Abell 2163.. .... .... 0,20 229±}f 15.6±f4 8.7+if 7.5- ! J 8.03±|^ 

Abell 2218.......... 0.18 206+{'( 10.0±}g ■ 4.8±g : g 4.1+g ; g 3.20±J ; |g 

Abell 2261 6,22 146±£? 6.4±j| 2.9±g ; | 2.5±g : | 0.76±g : |f 

Abell 267 0.23 .138+? 6.5+Jg 2.7±g : | 2.3±g | 0.75±g ; |? 

Abell 370 0.38 96+f 7.3±g : ? 3.0±g : | 2.6±g| 0.45±g (g 

Abell 665... ........ 0.18 181+^ 7.4+?? 3.3±g : | 2.8±g : g 1.55±g : |f 

Abell 697 0.28 130+Jg 8.3±? : § 3.7±g ; g 3.2±J : g 0.99±gi| 

Abell 773 0.22 150±§ 6.9+gg 3.1±g| 2.7±g'f 0.93±gif 

CL J0016+1609 0.54 81±§ 12.6±gg 4.1±g4 3.5±g ; | 0.92±g : JJ 

Cl, .) 1226-1-3332 0.89 53±? 9.2±gg 3.2±gf 2.7±g : l 0.27±g : gf 

MACS J 0647.7+7015 0.58 72±| 8.6±j^ 3.4±g : | 2.9±g : | 0.38±g : J| 

MACS J1311.0-0310. 0.49 73±| 6.8±gf 2.4±g ; f 2.1+gj 0.26±g : gf 

MACS J2214.9-1359. 0.48 91±| 10.2±g? 4.5+gf 3.9±gg 0.74±g : || 

MACS J2228.5+2036 0.41 100+f 11.0±} ; | 4.3±g : S 3.7±g : | 1.03±g : f? 

MS 0451.6-0305..... 0.55 84±g ■ 11.3±J : § 4.7±g : | 4.1±gf 0.79±g : ?£ 

MS 1137.5+6625.... 0.78 49±| 6.4±g| 1.9±g : f 1.6±g : | 0.12±g : g| 

MS 1358.4+6245 ... . 0.33 99±| 6.4±g : g 2.4±g j 2.1±g : f 0.40±g : g 

RX .11347.5-1145.... 0.45 108+4 12.7+JJ 6.5±g : ^ 5.6±g : g 1.41±g : f? 

RX J2129. 7+0005 ... . 0.24 123+^ 7.9±|? 2.1±?g 1.8±g : f 0.80±g : gf 

ZW 3146. — 0.29 128+g 8.3±£g 3.9+gf 3.3±gf 0.94±g ; l| 
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M tot( M ©) 

Fig. 2. — Scaling relations between the Sunyaev-ZeFdovich Effect Y and M gas M tot and kT e . 
Open squares (in black) are clusters at 0.14 < z < 0.30, open diamonds (in red) are clusters 
at 0.30 < z < 0.89. All measurements follow simple power-law models with indices that are 
consistent with the values of the self-similar scaling theory (Table 1). 
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Fig. 4.- — Comparison between the Sunyaev-Zel’dovich Effect Y parameter and the X-ray 
quantity Yx (see Section 5 for explanation of the normalization constant); open squares (in 
black) are clusters at 0.14 < z < 0.30, open diamonds (in red) are clusters at 0.30 < z < 0.89. 
The dashed black line corresponds to y = x. 
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Fig. 5. — M gas vs. Y for simulated and observed clusters. Open squares (in black) are 
OVRO /BIMA /Chandra measurements (as in Figure 2). In red (open and filled circles), 
simulated clusters from a cooling and star-formation feedback model; in green (open and 
' filled triangles) , simulated clusters from a non-radiative model (from Nagai 2006). Open 
symbols represent simulated clusters at z=0, filled symbols represent simulated clusters at 
z=0.6. 




